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ABSTRACT 

We present a spectroscopic study of the stellar cluster population of M82, the archetype starburst galaxy, 
based primarily on new Gemini-North multi-object spectroscopy of 49 star clusters. These observations con- 
stitute the largest to date spectroscopic dataset of extragalactic young clusters, giving virtually continuous 
coverage across the galaxy; we use these data to deduce information about the clusters as well as the M82 post- 
starburst disk and nuclear starburst environments. Spectroscopic age-dating places clusters in the nucleus and 
disk between (7, 15) and (30, 270) Myr, with distribution peaks at ~ 10 and ~ 140 Myr respectively. We find 
cluster radial velocities in the range vr G (—160, 220) kms" 1 (wrt the galaxy centre) and line of sight Nal D 
interstellar absorption line velocities u^ aID 6 (—75, 200) kms -1 , in many cases entirely decoupled from the 
clusters. As the disk cluster radial velocities lie on the flat part of the galaxy rotation curve, we conclude that 
they comprise a regularly orbiting system. Our observations suggest that the largest part of the population was 
created as a result of the close encounter with M8 1 ~ 220 Myr ago. Clusters in the nucleus are found in solid 
body rotation on the bar. The possible detection of WR features in their spectra indicates that cluster formation 
continues in the central starburst zone. We also report the potential discovery of two old populous clusters in 
the halo of M82, aged > 8 Gyr. Using these measurements and simple dynamical considerations, we derive a 
toy model for the invisible physical structure of the galaxy, and confirm the existence of two dominant spiral 
arms. 

Subject headings: galaxies: evolution — galaxies: individual (M82) — galaxies: kinematics and dynamics 
— galaxies: spectroscopy — galaxies: starburst — galaxies: star clusters 



1. INTRODUCTION 

The evolution of galaxies and their star formation rates are 
closely linked to interactions with other systems. The M81 
group of galaxies provides a good local laboratory for study- 
ing such interactions and the ensuing violent episodes of star 
formation. The irregular galaxy M82 is believed to have suf- 
fered a close encounter with M8 1 ~ 220 Myr ago (Brouillet et 
al. 1991; Yun, Ho & Lo 1994; Yun 1999). This interaction ev- 
idently triggered the powerful starburst event that is currently 
seen in the central region. The first systematic study of the 
centra l starburst was undertaken by lO'Connell & Manganol 
d 19781) who identified a number of optically high surface 
brightness starburst clumps (lettered A-E). These have sub- 
sequently each been found to contain hundreds of supe r star 
clusters (SSCs; |P ; Connell et al.|[l99l iMelo et al1l2005l) . Un- 
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under a cooperative agreement with the NSF on behalf of the Gemini partner- 
ship: the National Science Foundation (United States), the Science and Tech- 
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(Canada), CONICYT (Chile), the Australian Research Council (Australia), 
CNPq (Brazil) and CONICET (Argentina). 

- Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. These observations are associated with pro- 
gram #10853 

3 Department of Physics and Astronomy, University College London, 
Gower Street, London, WC1E 6BT, UK; isk@star.ucl.ac.uk 
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fortunately, however, the nearly edge-on inclination of the 
galaxy and its high dust content mean that, in projection, M82 
presents a highly irregular, confused mo rphology with optical 
obscuration levels reaching 10-15 mags (ISatvapal et al.l 19951: 
lAlonso-Herrero et al.l2003l) . 

Studies of star clusters outside of the central regions have 
been rel atively scarce because of t he hi gh levels of obscu- 
ration iGallagher & Smith! (1 19991) and ISmith & Gallagher! 
(1200 ll) studied the isolated star cluster M82-F using a combi- 
nation of HST imaging and ground-based spectroscopy. The 
derived age and mass of 60 ± 20 Myr and 1.2±0.1xl0 6 M Q 
(revisited in lBastian et al.ll2007h for M82-F show that signifi- 
cant star formation events have oc curred in the M82 di sk be- 
fore the current nuclear starburst. Ide Griis et al.l d2001l) stud- 
ied the cluster population of regi on B (a ~ 1 kpc regio n in the 
eastern part of the disk; see also Ide Grijs et al.ll2003l) . Using 
optical BVI photometry measured from HST/WFPC2 imag- 
ing, they found that the region B clusters have ages of 0.5-1.5 
Gyr, leading to region B being dubbed a 'fossil starburst' re- 
gion. However, not only are optical studies significantly af- 
fected by high levels of extinction, there also exists a degen- 
eracy between the integrated optical light properties of young 
(<300 Myr) and old (>1 Gyr) clusters when only using BVI 
photometry, since the baseline provided does not include the 
age-sensitive Balmer jump. This feature lies at ~ 3700 A 
and is therefore 'shared' between the U and B bands. It is 
thus difficult to break this degeneracy withou t imaging in at 
least one additional filter (lAnders et al.ll2004l) . This is due to 
the entanglement of the age and extinction measurements that 
have to be performed simultaneously when using photometry. 
In the absence of J7-band coverage, the simple stellar popu- 
lation (SSP) model tracks to which we compare overlap, giv- 
ing rise to a number of solutions - often vastly different with 
respect to age. The inclusion of U -band goes a long way to- 
wards breaking this degeneracy, producing a single solution, 
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i. e. a single pair o f age and extinction values. In our recent 
studie s of region B (ISmifh et al.ll2007b iKonstantopoulos et alJ 
120081 hereafter Paper I) we used HSTI ACS UBVI photometry, 
which covers the Balmer jump region, and Gemini-GMOS 
spectroscopy, where no degeneracy is present. Based on these 
data we showed that the region B clusters are much younger 
than previously reported, with ages peaking at ~ 150 Myr, 
and that region B is optically bright because it represents a 
window into the body of the galaxy, where the foreground ex- 
tinction is lower th an in the surrounding areas. A recent spec- 
troscopic study by iMavva et al.l d200 6) also found the stel- 
lar population of the M82 disk to be consistent with having 
formed in the last Gyr, during a single, short (~ 0.3 Gyr) 
burst of star formation. 

Cluster spectra are affected to a lesser extent by extinc- 
tion, and provide more reliable age determinations, as well 
as access to a wealth of other information (e.g. velocity 
and metallicity) — a lbeit at the cost of observing time (e.g. 
iTrancho et alj|2007b . In this paper we present optical spec- 
troscopy for 49 clusters in M82, comprising the largest spec- 
troscopic sample to date of young extragalactic star clusters. 
The observations were planned in such a way as to provide not 
only a statistically viable sample of the post-starburst cluster 
population, but almost continuous coverage across the disk 
and nucleus. This enables us to study the star clusters individ- 
ually and as a population, and also infer their parameters as a 
function of the M82 starburst environment. 

In this study we have sought to establish an age-dating 
technique based on the fit between observation and a range 
of models, that is as free as possible from the effects of ex- 
tinction (which, as alluded to above, is a significant obstacle 
in M82). To diagnose the cluster properties, the main tool 
we employ is the optimised fitting of evolutionary synthesis 
models (in this case SSPs) to both the optical cluster spec- 
tral energy distributions (SEDs), and to individual Balmer ab- 
sorption lines. This automatically breaks the aforementioned 
degeneracy between young/old ages, since the shape of the 
Balmer absorption lines is very sensitive to age. Furthermore, 
effects due to emission from foreground gas can be mitigated 
by careful masking of the line profile. 

This paper is organised as follows: we present the imaging 
and spectroscopic datasets in § 2; the measurements of ages 
and radial velocities from spectroscopy, and colours and red- 
dening from photometry are presented in § 3; we discuss the 
implications of these results in § 4; and propose a model for 
the obscured physical structure of M82 in § 5. We summarise 
all our results in the final section, § 6. We adopt a distance to 
M82 of 3.6 Mpc (Freedman et al. 1994). 

2. DATA ACQUISITION AND REDUCTION 
2.1. The full M82 spectroscopic dataset 

The results presented in this paper are based on spec- 
troscopy obtained with the Gemini-North Multi-Object Spec- 
trograph (GMOS-N) on the Gemini-North telescope, as part 
of observing program GN-2006A-Q38 and the follow-up pro- 
gram GN-2007A-Q21 (PI L. J. Smith; hereafter Q38 and 
Q21). Overall, two multi-object slit masks were designed and 
included a total of 54 slits, from which we extracted 62 spectra 
in the range AA3700-6500 A. Figure Q] marks the positions of 
the slits on HSTI 'ACS (Advanced Camera for Surveys) imag- 
ing of M82 dMutchler et al.ll2007l) and presents all the cluster 
spectra used in the subsequent analysis. 

The source selection was performed within the bounds of 



the original observing program (Q38). The candidate young 
massive clusters (YMCs) were selected based on GMOS pre- 
imaging after cross-identification with HST imaging, with the 
aim of choosing the brightest isolated clusters, spanning most 
of the optical extent of the disk. Starting from a masterlist 
of several hundred candidates, we built two multi-slit masks 
(one for each observing program), with a standard 0'.'75 slit 
width, while the slit length was varied to optimize the overall 
arrangement. In all cases, we attempted to apply a minimum 
slit length of 6'.' 5, in order to provide adequate sky for back- 
ground subtraction. In Q38, we fell below 6'.'5 in a few cases, 
to accommodate more clusters along the slit; this, however, 
complicated the aperture extraction procedure. For that rea- 
son, in Q21 we did not sacrifice the slit length in favour of 
cluster number, particularly since the mean magnitude was 
fainter compared to Q38. 

The GMOS observational setup was identical for Q38 and 
Q21. We used the B600 grating (0.9 Apx -1 ) and the three 
GMOS CCDs in 2 x 2 binning mode. The observations were 
executed in pairs of pointings with different central wave- 
lengths (5080 A and 5120 A) in order to cover the inter-chip 
gaps. A brief journal of observations is given in Table[T] The 
total exposure time for each dataset was eight 1800 s expo- 
sures for a total on-source integration time of 4 hr. CuAr arc 
exposures and Quartz-Halogen flat fields were taken in be- 
tween target exposures, and bias frames were taken as part of 
the Gemini base-line calibrations (GCAL). 

The data reduction was performed using standard IRAF 9 
and purpose-designed Gemini-IRAF routines. In brief, this 
consisted of bias-subtraction, flat-fielding, and the combina- 
tion of a three-CCD image into a single frame, while cor- 
recting for the variation in detector quantum efficiency. We 
used the obtained CuAr frames for the wavelength calibra- 
tion, which resulted in satisfactorily low residuals (typically 
less than 0.2 A). As Gemini does not have a functioning 
atmospheric diffraction corrector, we corrected for differen- 
tial atmospheric refraction using a purpose-built IDL (Inter- 
active Data Language) routine written by B. W. Miller (pri- 
vate communication; this routine is based on lFilippenkol 19821 
where the effects of atmospheric differential refraction are 
discussed). Finally, all eight available exposures of each tar- 
get were extracted and then flux calibrated using exposures of 
standard stars Wolf 1346 (observed 2006 May 25) and HZ 44 
(observed 2008 Feb 18) for Q38 and Q21 respectively. 

In summary, the two masks contained 54 slits (see Table [B 
and the total number of sources observed was 62, of which 28 
have a S/N above 20 (conservatively measured in two contin- 
uum windows covering AA4500-4800 A and AA5000-5500 
A). We present spectra for 49 clusters in this paper, and omit 
the remaining 1 3 on grounds of low data quality. The spectral 
resolution of the Q38 and Q21 cluster sample was measured 
to be 3.5 A and 3.7 A respectively from the CuAr spectra, but 
we assume a common resolution of 3.7 A for the sake of sim- 
plicity. The spectral range of our data is « AA3700-6500 A, 
which includes the H/3 and H7 lines, our main age indicators, 
and the Na I D doublet, which we use as a tracer of interstel- 
lar gas radial velocities. The exact wavelength range varies 
with each slit, depending on its positioning with respect to the 
M82 major axis, which defines the dispersion axis. Our ef- 

9 IRAF is distributed by the National Optical Astronomical Observato- 
ries, which are operated by the Association of Universities for Research in 
Astronomy, Inc. under contract with the National Science Foundation 
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FIG. 1 . — The top panel shows the existing dataset of star cluster spectroscopy (programmes Q38 and Q21, 16 hours of GMOS-MOS on-target integration time 
in all), with slit positions plotted over HST-ACS imaging of M82 (Hubble Heritage ACS mosaic of 8 WFC images, Mutchler et al. 2007). The dataset samples 
individual star clusters across the entire galaxy disk and cluster associations in the nucleus. In the bottom panel we present all spectra used in this study, grouped 
according to the region the clusters occupy: 'W, 'B', 'N' and 'E' stand for western disk, region B, nuclear region and eastern disk respectively. Note that we 
have excluded spectral regions with S/N < 20 from the plot. 



fective wavelength range for all slits is AA4000-6000 A; this 
is set by the poor sensitivity of GMOS below 4000 A and 
an uncertain flux calibration above 6000 A. Fig. |2] demon- 
strates these features on the representative high S/N spectrum 
of cluster 112.2. Positional and other basic information for 
all clusters can be found in Table [2] (the coordinates of our 
sources have been converted to the HST WCS, as the mosaic 
imaging is publicly available). 
The cataloging convention we have adopted for cluster 



naming follows the increasing slit number in the masterlist 
compiled during the initial source selection. The list is not 
continuous (as more sources were initially selected than could 
be accommodated), but the increasing number does indicate 
the position of a cluster along the M82 disk - numbers count 
from west to east, along the major axis. In addition, several 
slits include multiple sources; in these cases, we refer to clus- 
ters by the slit number with a decimal part, to indicate the 
increasing number of the source (again, counting from west 
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FIG. 2. — A representative flux-calibrated, high S/N spectrum, with the 
main absorption features denoted. 



TABLE 1 

Summary of spectroscopic observations taken with 
Gemini-North 



Dataset 


Dates of observation 


Seeing" 


No. slits 


Exposure time 


GN-2006A-Q38 


2006 April 5 


O'.'S 


39 


8 x 1800 s 


GN-2007A-Q21 


2007 June 4-8 


0V45 b 


15 


8 x 1800 s 



to east). Our sample spans most of the extent of the optically 
visible disk, which is approximately 6 by 2 kpc, as seen in 
HST imaging. 

2.2. HST Photometry 

We also make use of imaging data of M82 to derive clus- 
ter colours and deduce information about their surroundings. 
We use the publicly available Hubble Heritage ACS mosaic 
imaging of M82 (for a full description, we refer the reader to 
iMutchler et al.ll2007l) . which covers the galaxy disk, nucleus 
and outskirts in six pointings and in three filters, F435W, 
F555W, F814W. These are roughly equivalent to the stan- 
dard Johnson BVI, so we use this notation throughout this 
paper. Note, however, that we do not perform any conver- 
sion between the two systems. All of our observed clusters lie 
within the field of view of the ACS imaging. The spatial scale 
of the images is 0.05 arcsec px _1 , and one pixel corresponds 
to ~ 0.9 pc at the distance of M82. 

We used the HST images to perform aperture photometry 
on all sources in our spectroscopic sample. For most disk 
clusters, we applied an aperture radius of 10 pixels to avoid 
contamination by nearby objects, with a 2 pixel wide sky an- 
nulus placed at 12 pixels. As most of our sources are of con- 
siderable size, we calculated aperture corrections for each fil- 
ter as —0.34 mag for the B and V bands and —0.45 mag for 
the / band, based on the flux difference between 10 and 30 
pixel apertures for isolated clusters. After close inspection of 
each source individually, it was deemed necessary to limit the 
aperture size in some cases, due to the presence of contami- 
nating neighbouring sources. Thus, we applied 5 pixel aper- 
tures with the same annulus for eight disk clusters as well as 
all sources in the nuclear region (in the nucleus, the sky annu- 
lus was placed between 7 and 9 px), and calculated the correc- 
tions to 30 pixels as -0.94, -0.90 and 1.04 mag for the B, V, 
and I bands respectively. While these aperture corrections are 
significant and dependent on cluster size, the cluster colours 
should be largely unaffected. The photometric results are pre- 
sented in Tab le [H this includes cluster reddening, which is 
discussed in § 13.21 



3. SPECTROSCOPIC AND PHOTOMETRIC ANALYSIS 



3.1. Cluster ages and radial velocities from spectroscopy 

In this section we present our methodology for deriving 
cluster ages and illustrate our techinques through example 
fits. We employ the Cumulative x 2 Minimisation method 
(CCM) presented in Paper I to obtain ages for the full sam- 
ple of 49 star clusters in the disk and nuclear region of M82. 
This method is based on the goodness of fit between ob- 
served cluster spectra and evolutionary population synthe- 
sis (EPS) models. We chose t o compare our spectra to the 
iGonzalez-Delgado et al.l d2005L hereafter GD05) model li- 
brary, a range of theoretical models specifically compiled 
for young star clusters. The major advantages these models 
present over more establ ished empirical cluster SEDs (e. g. 
iBruzual & Charlo3l2003l) are the high spectral and temporal 
resolution (0.3 A and a few Myr for the first 100 Myr). We 
used the full range of models (4 Myr to 17 Gyr) compiled us- 
ing Padova isochrones, and a Salpeter IMF. We have assumed 
a sola r metallicity following the findings of iMcLeod et alj 
(j!993j) for the M82 ISM. 

Before the fit can take place, we rectify both observed spec- 
trum and template (using two continuum windows adjacent 
to the absorption line) and smooth the models by a factor of 
three and then re-bin them to the wavelength scale of the ob- 
served spectra. In order to obtain a first estimate of the cluster 
age, we apply an approximate radial velocity (RV) correction 
on the observed spectra (matching the Balmer line profiles by 
eye; this calculation is later refined). We use the lower Balmer 
series line profiles (H/3 and H7) to perform the fits, and take 
care not to include any emission from ionised gas in the line 
of sight (which is plentiful in this gas rich, edge-on galaxy). 
This makes our age determination technique largely indepen- 
dent of the effects of foreground gas and dust. 

In brief, the CCM calculates the \ 2 f° r eacri point in the 
fit, and picks the best fitting age as the model with the low- 
est 'reduced' x 2 , expressed through x 2 = ^X 2 /d (note that 
we will omit the subscript v hereafter); d is the number of 
degrees of freedom, taken as the number of points included 
in the fit. This appr oach follows the reasoning outlined in 
lLampton et al.1 d 19761) . Our fitting technique is therefore sen- 
sitive to the overall line profile (core and wings) rather than 
just the line depth. This approach overcomes both the line 
strength degeneracy (discussed in § [U that exists between 
models of age - 200± 100 Myr and - 600± 100 Myr, and the 
lack of information on the line core where the superimposed 
emission is strong or complex. We therefore obtain two mea- 
surements of the cluster age, based on the fits on the H7 and 
H/3 features. As the two measurements do not always agree 
perfectly (usually due to the lower S/N in the H7 spectral re- 
gion), we examine each case individually. In a few cases, one 
of the two measurements is excluded from the fit; we apply 
this principle by enforcing a minimum acceptable S/N ratio 
of 10 (this applies to 12 clusters in the sample). 

The uncertainty on a given fit is defined as the range of 
ages that lie within la of the minimum \ 2 - quantified by 
the criterion yJ < Y^Lv, + N, where N is t he number of 
free parameters in the fit ( lLampton et aljfl976l) . In this case, 
the only varying parameter is the age of the model, therefore 
N = 1. This process is slightly complicated by the fact that 
a large number of our adopted ages consist of the statistical 
mean of the results of two x 2 fits. In these cases, we combine 
the probability distributions derived for each fit (simply by 
adding them) and use a value of N = 2 to define the la confi- 
dence region. We stress that, while this error determination is 
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TABLE 2 

Positional information (equinox J2000) and aperture photometry for all 49 clusters in the spectroscopic sample. 



ID (alt ID) 




RA 






DEC 




TUB 


± 


mv 


± 


mi 


± 


A v 


M v 


± 




(h 


m 


s) 


(° 


/ 


") 


(maj 


;) 


(mag) 


(mag) 


(mag) 


(mag) 


7 


9 


55 


27.53 


69 


40 


3.7 


21.02 


0.09 


20.55 


0.07 


19.67 


0.05 


0.85 


-8.09 


0.95 


8 


9 


55 


29.99 


69 


39 


53.6 


21.22 


0.10 


20.70 


0.08 


19.79 


0.06 


0.88 


-7.96 


0.91 


14 


9 


55 


35.46 


69 


39 


56.5 


21.93 


0.13 


20.89 


0.08 


19.48 


0.05 


2.03 


-8.92 


0.44 


17 


9 


55 


35.28 


69 


40 


18.5 


20.70 


0.08 


20.08 


0.06 


19.08 


0.05 


1.10 


-8.80 


0.80 


18 


9 


55 


35.99 


69 


40 


14.9 


20.37 


0.07 


19.76 


0.06 


18.74 


0.04 


1.15 


-9.17 


0.80 


19 


9 


55 


37.11 


69 


40 


20.9 


20.76 


0.08 


20.24 


0.06 


19.39 


0.06 


0.98 


-8.51 


0.87 


20 


9 


55 


37.50 


69 


40 


26.0 


21.25 


0.10 


20.77 


0.08 


19.94 


0.06 


0.75 


-7.76 


1 .03 


21 


9 


55 


38.56 


69 


40 


27.0 


21.01 


0.09 


20.31 


0.07 


19.41 


0.06 


1.02 


-8.49 


0.83 


28 


9 


55 


41.03 


69 


40 


27.8 


21.27 


0.10 


20.46 


0.07 


19.16 


0.04 


1.70 


-9.02 


0.53 


29 


9 


55 


41.85 


69 


40 


22.1 


21.07 


0.10 


20.07 


0.06 


18.55 


0.04 


2.17 


-9.88 


0.46 


33 


9 


55 


42.40 


69 


40 


37.8 


21.54 


0.11 


20.26 


0.06 


18.46 


0.03 


2.80 


-10.32 


0.37 


34 


9 


55 


43.69 


69 


40 


23.4 


21.11 


0.09 


19.82 


0.06 


18.14 


0.04 


2.62 


-10.58 


0.40 


38 


9 


55 


45.12 


69 


40 


36.6 


20.97 


0.09 


20.20 


0.06 


19.09 


0.05 


1.42 


-9.00 


0.63 


43.1 


9 


55 


45.61 


69 


40 


40.1 


20.63 


0.08 


19.14 


0.05 


17.13 


0.02 


3.22 


-11.87 


0.37 


43.2 


9 


55 


46.62 


69 


40 


43.9 


19.90 


0.06 


19.27 


0.05 


18.35 


0.05 


1.02 


-9.54 


0.93 


51 (F a ) 


9 


55 


47.08 


69 


40 


42.3 


17.61 


0.02 


16.50 


0.01 


14.94 


0.01 


2.42 


-13.70 


0.57 


58 


9 


55 


50.27 


69 


40 


35.0 


20.25 


0.08 


19.72 


0.07 


18.75 


0.07 


1.12 


-9.19 


0.82 


67.1 


9 


55 


50.44 


69 


40 


47.2 


19.92 


0.05 


18.72 


0.04 


17.04 


0.05 


2.47 


-11.54 


0.47 


67.2 


9 


55 


50.77 


69 


40 


47.9 


21.18 


0.13 


20.69 


0.26 


18.78 


0.15 


2.17 


-9.27 


0.44 


67.4 


9 


55 


51.93 


69 


40 


50.2 


21.60 


0.12 


20.80 


0.1 1 


19.12 


0.22 


2.12 


-9.11 


0.43 


78.1 


9 


55 


52.64 


69 


40 


46.9 


21.16 


0.14 


20.76 


0.23 


19.50 


0.40 


1 .20 


-8.22 


0.69 


78.2 


9 


55 


53.13 


69 


40 


49.1 


20.43 


0.10 


19.36 


0.08 


16.89 


0.07 


3.45 


-11.87 


0.35 


91 (H b ) 


9 


55 


54.61 


69 


41 


1.6 


19.20 


0.04 


18.03 


0.02 


16.44 


0.01 


2.40 


-12.15 


0.51 


97 


9 


55 


55.65 


69 


41 


5.9 


19.08 


0.04 


18.56 


0.03 


17.63 


0.02 


0.90 


-10.12 


1. 12 


98 


9 


55 


55.50 


69 


41 


9.9 


19.69 


0.05 


19.31 


0.04 


18.58 
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not statistically flawless, it provides an adequate description 
of the range of models that can plausibly describe the spec- 
tra. This range therefore matches what one might pick 'by 
eye'. This statistica l analysis was also successfully used in 
iBastian et all (120081) . 

Our fitting method is demonstrated in Fig.[3]for three clus- 
ters of varying luminosity and S/N. We present a BVI colour- 
composite image (the spatial scale is a square of side 200 px 
or 175 pc) and the spectrum of each cluster, with the best- 
fit model overplotted. We also show the probability distri- 
butions of all available models as a function of age below 
each spectrum. We note that, although not demonstrated, the 
H7 and H/3 fits do provide overlapping \ 2 distributions (i. e. 
equivalent fits). In the top panel we show a low S/N spec- 
trum (cluster 19, S/N < 20), where the H/3 region provides 
the only reliable fit. The middle panel shows the H7 profile 
in a medium S/N (~ 25) spectrum. For this cluster, we are 
able to obtain two age measurements in good agreement. In 
the bottom panel we show the high S/N (~ 40) spectrum of 



cluster 112.2, where the H/3 and H7 measurements give the 
same age. We note that these plots do not present a subset of 
high-quality fits, but are meant to provide an overview of all 
problems faced by our fitting method, be they due to low S/N, 
embedded/blended emission, or related to the template fitting 
process. 

One immediately noticeable feature in these fits is the ex- 
istence of two best-fit troughs. While the two fits are statisti- 
cally equivalent, there is a significant volume of evidence that 
favours the younger ages. This is based on the comparison 
to ages independently derived through UBVI photometry, as 
well as a number of considerations of the star formation his- 
tory of M82 that will be discussed in § 14.31 The age measure- 
ments are shown in Table [3] where we have listed the accept- 
able range as the extent of the preferred best fitting trough. 

A final consideration concerns the effect of our choice of 
model library. We used the GD05 templates to achieve the 
highest quality fits, given their high spectral and temporal res- 
olution. As a consistency check, we fitted some of the high 
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FIG. 3. — (Previous page). Cluster age dating. Each of these three image-spectrum pairs represents a BVI composite image and part of its optical spectrum, 
as used in our age-dating method. In the images, East is to the left and North towards the top, with the major axis of the galaxy inclined at ~ 30 degrees to 
the horizontal. The spatial scale is a square size of 200 px, or 175 pc. In the spectra, the top panel shows the fit, with the solid line and dashed (green) lines 
representing the observed spectrum and best fitting model respectively. The dashed blue boxes indicate the spectral regions where the fit takes place. The bottom 
panel shows the probability distribution of the CCM across age-space, with a vertical line indicating the best fitting age, while the horizontal line indicates the 
error on the fit (Xmin + x )- 

Top: Cluster 19, in the western side of the galaxy disk. This is an example of a low S/N spectrum (S/N < 20 at 5000A). The faint and strong absorbing dust 
lanes to the west and east of the cluster provide a demonstration of the sub-parsec spatial scale and the amount and strong effect of the galactic dust in M82. This 
cluster has two high probability regions in age-space, however, the CCM draws its best fit from the low age trough, placing the age at ~ 40 Myr. 
Middle: Medium S/N (~ 25) spectrum of cluster 38. Note that the fit presented here is on the H7 line, and the S/N is lower at the blue end of each spectrum. 
Bottom: High S/N (~ 40) spectrum of cluster 1 12.2, situated north of region B. Note that in this case, all age measurements (for both lines) agree on 199 Myr. 
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S/N c lusters in our sample with BC03 dBruzual & Charloll 
120031) models of similar characteristics (same isochrones, 
metallicity, and IMF) and obtained very similar results. For 
instance, cluster 91 was found to have an identical age under 
GD05 and BC03. Further in support of our choice of tem- 
plates is the fact that, under BC03, the best fit troughs are 
represented by very few consecutive models (perhaps three 
or four), whereas with GD05 the x 2 probability distribution 
appears smooth and therefore more reliable. 

Having measured the age of each cluster, we can now de- 
termine their radial velocities (RVs). We used the IRAF FX- 
COR routine to cross-correlate the best fitting model with 
the cluster spectrum. We also manually measured the ve- 
locities of the Na I D lines (AA5890, 5896) which originate 
from cold interstellar gas along the line of sight (LOS). We 
transformed the measured cluster velocities to the heliocentric 
frame of reference and corrected for th e systemic velocity of 
M82 (200kms- 1 : lMcKeith et alHl993h . The final RV values 
vary between — 220km s~\ with associated Nal D veloc- 
ities in the range 5 — 200 km s -1 (in absolute value). These 
measurements are also presented in Table [3] We also used 
the radial velocities to refine our age measurements: we re- 
peated the age fitting process, this time applying a more accu- 
rate Doppler-shift according to the derived RVs. This allowed 
for a better alignment between spectrum and model, which 
had been estimated in the first iteration. This second iteration 
of the age fits did not alter any of the best fitting ages. 

We divide our sample into disk clusters and nuclear region 
knots, where the nuclear region is defined as a projected ra- 
dius of_£^_2J)0_2cabj3ut the dynamical centre of the galaxy 
(see iMcKeith et al1l!993l for exact location). Note here that 
we use the terms 'cluster', 'knot' and 'cluster complex' inter- 
changeably for these nuclear region objects; they are proba- 
bly not individual clusters but cluster complexes (physically 
associated) or asterisms (chance projections). We find all disk 
clusters to have ages between ~ 30 and 270 Myr, with a mean 
(and median) value of ~ 140 Myr. For the knots in the nuclear 
region (slit numbers 67 and 78), these values range from 7 to 
30 Myr with a mean of 10 Myr. Unfortunately, it is not pos- 
sible to obtain accurate age dating in the nuclear region as the 
fits are dominated by strong emission lines. This restricts the 
fitting to a spectral range that hardly covers the Balmer ab- 
sorption lines. The CCM picks very young ages due to the 
flatness of these spectra; in the case of knots in slit 67, the 
fit favours the very youngest models in GD05 library, as all 
models above a few 10 Myr feature noticeably stronger metal 
lines 10 . Further evidence of the young age of these clusters is 
offered by the strong, embedded emission lines. Such strong 
features centred at the same wavelength as the absorption can 
only be associated with the clusters. The clusters in slit 78 dis- 
play a similar appearance, however, with some more structure 
(i. e. they are not as 'flat' as the sources in slit 67) thus they 
are assigned slightly older ages (in the third and fourth decade 
of Myr). In addition, we find offset emission in cluster 78.1 
perhaps indicating that it is not situated deep in the galaxy like 
the rest of the 'nuclear' knots. Finally, our sample also con- 
tains two halo objects, which we have dated at > 8 Gy r; we 
discuss these candidate old populous clusters (OPC) in § 13.31 

3.2. Cluster photometry: colours and extinction 

10 The fit extends to the regions adjacent to the Balmer lines, where weak 
metal lines are present. 



2.0 




vijctors: 












- 




_ 1 ; S My 
_ 2; 16 My 


3mag 
4mag 


f 


57 

1 






1 135 

f 




1.5 


3: 248 My 

_ 4: 3.2 Gyr 
. 5: 6.0 Gyr 


4mag 
0.8 mag 
2mag 






56 




3 

2 . 




















E 

^ 1 -° 








6.0 C 
3.2 Gy 








_ZB.2 


0.5 

0.0 


* < Vk-r 


550 Myr 


^12 My 




\ 


_135 







0.0 0.5 10 1.5 2.0 2.5 3.0 

V — I / mag 



FIG. 4. — B — V vs V — I colour-colour diagram for all disk clusters with 
available spectroscopy; we have marked off-disk OPC candidates with stars 
and nuclear knots with triangles. Foreground-corrected cluster extinctions 
range between near zero and 2.7 magnitudes (Ay) in the disk and 1.4 — 
3.0 mag in the nuclear region. The concentration of datapoints around the 
250 Myr extinction vector reflects the limited age scatter about that age. The 
two indicated pairs of isolated clusters (nos. 113, 135, 156 and 157) have 
underestimated fluxes in one or more of the three bands. 

We use the photometric measurements to derive clus- 
ter colours (B — V and V — I) and present these re- 
sults in Fig. [4] Here we compare the observed clus- 
ter colours with GALEV evolutiona ry models for SSPs 
( I Anders & Fritze-v. Alvenslebenl 120031) . As for our chosen 
SEDs, the photometric models use the Padova isochrones, a 
Salpeter IMF and solar metallicity. 

This plot can be used as a rough diagnostic of age but by 
no means an accurate one, given the lack of information in 
the Balmer jump region and the consequent degeneracy in the 
BVI baseline. From this it becomes obvious that there is no 
clear distinction when age-dating photometrically between a 
cluster of, e. g. age 4 Myr and 2.5 Gyr. Given some a priori 
knowledge of extinction, only some rough deductions can be 
made. 

More importantly, given a strong prior knowledge of age 
(which we possess from spectroscopy), colours can be used 
to estimate the reddening for each cluster. To do this, we used 
a simple code to transfer each cluster colour along the extinc- 
tion vector until it is close to the intrinsic (model) colour at 
the derived spectroscopic age. Assuming a standard Gal actic 
extinction law with Ry =3.1 dRieke & Lebo fskv 1985|), we 
measured extinctions in the M82 galaxy disk in the Ay range 
0.4-3.2, with a median value of 1.3 mag 11 . The correspond- 
ing range of values for the nuclear region is 1.4-3.5 with a 
median of 2.2 mag. To derive approximate errors for the Ay 
values, we have adopted an age range of 50 Myr (a represen- 
tative value for the error in age measurement) and calculated 
the corresponding extinctions. 

All values are listed in Table 12 Note that the foreground 
extinction from the Milky Way is ~ 0.5 mag in the F-band, 
so the lowest values we derive correspond to roughly no ex- 
tinction from material associated with M82. Interestingly, the 
measured values do not reflect the large amounts of extinction 
one might expect to find in this galaxy - reaching 15 magni- 
tudes, as we mention in the introduction. This is due to our 
sample, which sought out the brightest possible sources to se- 

" This range excludes cluster 158, a faint cluster with highly uncertain 
photometry 
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cure the acquisition of high S/N spectra. As a result, no highly 
extinguished sources were included. 

3.3. Observations of off-disk clusters 

Our spectroscopic dataset includes two halo objects located 
at a projected distance of ~ 400 pc to the south of the east- 
ern galaxy disk, selected from the initial masterlist as candi- 
date old populous clusters. We fit those objects (nos. 152 and 
139) with GD05 models of |Zq (a value characteristic of old 
globular clusters) and find best fit ages of > 8 Gyr. In the 
case of 152, we also find a peculiar/halo radial velocity as it 
lags significantly behind the disk population. Unfortunately, 
the low S/N of cluster 139 does not allow for a certain age fit 
or a high-confidence radial velocity cross-c orre lation, there- 
fore only 152 was investigated further (see § |4.61 l. These clus- 
ters have not been observed in the past (they are not included 
in the list by ISaito et al.l|2005l who studied a sample of GC 
candidates around M82 and reviewed the case for GCs in the 
galaxy) and make very interesting targets for further study. 

3.4. Detection of WR features in the nuclear region 

The designed multi-object mask features two long slits 
placed over the M82 nuclear region, covering a total of seven 
knots. The precise location of these sources is not easy to 
pinpoint, given the apparent complexity of the nuclear region. 
We therefore reserve the possibility that they lie near the sur- 
face of the galaxy and are p rojec ted onto the nuclear region 
(this matter is revisited in § |4.51 >- In addition, the extracted 
apertures are likely to contain blended light from more than 
one star cluster, therefore we do not treat these as observations 
of individual sources. The activity in the nuclear region is key 
to understandin g the st arburst history of M82. A recent study 
by iBarker et al.l d2008l) . found no evidence of cluster forma- 
tion over the last 4 — 6 Myr in the nuclear region. However, we 
find Wolf-Rayet (WR) features in the spectra of three clusters 
in the region (maybe cluster complexes; sources 1, 2 and per- 
haps also 4 within slit 67), which are indicative of such young 
ages. More specifically, we find evidence of He II A4686, C III 
A4651, and possibly Civ A4659, Nm A4634 - 4640 and 
N V A4603 - 4619; all these features form part of the 'blue 
WR b ump' at 4600 - 4700 A dSidoli et al.ll2006t iBibbv etaT] 
120081) . In Fig. [5] we provide a finding chart for slit No. 67 
and we present the spectra of apertures 1 and 2, indicating the 
locations of the proposed WR features. This detection needs 
to be confirmed with higher S/N in this spectral region, and if 
proven to be true, it will classify M82 as a WR galaxy. 

We stress at this point that the exact source of these emis- 
sion lines is uncertain. They could originate from (or be en- 
hanced by) scattered light from the nuclear region, rather than 
the observed knots. In order to avoid contamination by nebu- 
lar features in our analysis, we have extracted these spectra in 
narrow apertures. Despite having taken that precaution, we re- 
serve the possibility that the feature at A4658 is in fact nebular 
[Fe III]. In essence it is not the source of this emission that we 
report here, but simply its existence. Even though our data are 
not suitable for an in-depth analysis of the detected lines, they 
provide evidence that the nuclear starburst is still active. WR 
features only appear in the youngest stars (2 — 6 Myr), and are 
a reliable age-dating tool for young star clusters dSidoli et alj 
l2006t iBastian et aLl l2006). Mo re evidence for an ongoing nu- 
clear starburst was offered bv lSmith et al.l d2006l) . who found 
an age of 6.4 Myr for cluster Al in the nuclear region. The 
high mass of this cluster (~ 10 6 M ) implies the existence of 



a fully sampled cluster mass function at that era, and there- 
fore the contemporary formation of a large number of clusters 
- such a massive cluster represents but the 'tip of the iceberg' 
in terms of the underlying distribution. Hence we confirm that 
star/cluster formation is currently continuing in the nuclear re- 
gion ofM82. 

Additionally, it was proposed by IBarker et alj ( 120081) that, 
while stars are forming in the nuclear region of M82, star 
and cluster formation have become dissociated in this region. 
Our detectio n of young clusters in the same region studied by 
IBarker et al.l indicates that star and c luster formation is still 
ongoing. Furthermore. iBastianl d2008l) showed that cluster Al 
fits the relation between the brightest cluster in a population 
and the SFR of the host galaxy (a relation which extends from 
quiescent to extreme starburst galaxies). This suggests that 
star and cluster formation is proceeding normally in the nu- 
clear region - subject to our interpretation that these clusters 
reside in the nucleus (see £j |4.5l for further discussion). 

4. DISCUSSION 

In this section we discuss the spectroscopic (age and radial 
velocity) and photometric (colours and reddening) measure- 
ments for our sample of 49 disk, halo and nuclear star clusters 
in M82, and the implications of these results. 

4.1. Colours and reddening 

Cluster colours are presented in Fig.|H this plot works well 
in confirming the timescale for the starburst, with many dat- 
apoints distributed closely about the 250 Myr extinction vec- 
tor. In addition, the nuclear region clusters are intercepted by 
younger age vectors, and the OPC candidates by vectors in 
accord with their spectroscopic ages. There are two isolated 
pairs of clusters (nos. 113 and 135, 156 and 157) that do not 
agree with the spectroscopic ages; of those, cluster 113 ap- 
pears to have an underestimated flux in blue light, which we 
interpret as the result of differential extinction across the face 
of the cluster (this will shift the datapoint r edwards in blue 
bands preferentially; see Bastian et al. 2007, and Paper I for 
more discussion of such cases). Cluster 135 appears diffuse in 
the ACS imaging, and could therefore have disproportionately 
underestimated fluxes. Cluster 156 has an underestimated B- 
band flux, as it resides behind a thin, filamentary dust lane; 
cluster 157 lies near the edge of the optically visible disk, at 
the boundary between the bright, blue disk and wide, dark 
band to the north of region B, and is therefore differentially 
extinguished. 

We present cluster reddenings in the bottom panel of Fig.[HJ 
where in this case the values have be en corrected for fore - 
ground galactic extinction (Ay = 0.53; Sch legel et al.ll 1 998b . 
The values presented elsewhere in the paper do not incorpo- 
rate such a correction. 

4.2. Cluster Masses 

Having obtained accurate photometry (magnitudes and ex- 
tinctions) and ages for each star cluster in our spectroscopic 
sample, we derived their photometric masses (see Table[3]for 
results). We find average and median values of 2.5 x 10 5 and 
7.2 x 10 4 Mq. According to our photometric mass estimates, 
the most massive individual clusters in our sample are M82-F 
and H (nos. 51 and 91 on our list) and no. 152 (the candidate 
OPC), with values of 4.5 x 10 6 , 2.1 x 10 6 and 2.1 x 10 6 M . 
The total baryonic mass contained in the observed clusters 
sums to 1.1 x 10 7 M Q . 
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FIG. 5. — GMOS g-band imaging of the nuclear region of M82. Left: Positioning of spectrum extraction apertures in slit no. 67. The red box represents 
the exact positioning of slit 67 (tilted by 4.1° with respect to the major axis). We have overplotted the integrated light contained within the slit to identify the 
extracted objects, indicated by the dotted blue circles. Right: Extracted spectra for apertures 1 (black) and 2 (green/grey) in slit 67 showing the features, that 
usually form part of the 'blue WR bump' (~ 4600 — 4700A). We have normalised the spectra using two continuum windows just outside the plotted range, to 
show the WR bump. Note that, although we extracted these spectra using narrow apertures to avoid nebular contamination, we cannot reject the possibility that 
the line marked ' X ' is in fact nebular [Fe III] (as opposed to the C IV WR feature). 
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FIG. 6. — Cluster masses plotted against age (the age-mass diagram) for 
the entire spectroscopic sample presented in the current work. The solid line 
represents the detection limit at my = 21 (corresponding to My ~ —7 
at the distance of M82) translated into mass. The dynamical nature of this 
limit is due to the evolutionary fading of star clusters and therefore the vari- 
ation of the mass to light ratio (the line has been derived using GALEV 
models I Anders et al. 2004)). We have also shown the result of one mag- 
nitude of extinction as a dashed line (the same intrinsic limit coupled with 
one magnitude of extinction). The limit seems to accommodate all the data- 
points comfortably within the errors (not shown; of the order of 0.3 dex both 
horizontally and vertically); the dotted horizontal line shows a mass limit at 
~ 6 X 10 4 Mq, which is the mass to which our sample is predicted to be 
complete to an age of 1 Gyr; this is expressed as the point of intersection be- 
tween the mass cut and the detection limit. As a large fraction of our sample 
lies above that cut-off, we deduce that the non-detection of clusters older than 
~ 300 Myr is does not result from incompleteness), but reflects the intrinsic 
age distribution of the M82 starburst clusters. 



Based on these mass measurements, we can derive an esti- 
mate for the total mass in clusters for the entire galaxy. We 
generated an array of simple stochastic populations of star 
clusters, governed by a power law mass function of index —2, 
with a minimum mass of 50 M Q and no upper mass limit. As 
we know of five clusters in M82 with mass above 10 6 M©, we 
then varied the number of clusters in the population, N, un- 
til the average (over 100 runs) number of clusters with mass 
greater than 10 6 M Q is 5. From this we find a total cluster 
mass for M82 of ~ 10 8 M Q , the mass of a stochastic cluster 



population with N = 10 5 . This is an approximate lower limit, 
as there may exist more star clusters in M82 with mass greater 
than 10 6 Mq that have not been observed due to extinction. 

We repeated this exercise analytically, by calculating the 
missing mass according to a power law mass function of in- 
dex —2. Given five clusters with mass of the order of 10 6 , 
we integrate the mass function and extrapolate a mass of 
~ 6 x 10 7 M Q for the observable population. Taking this 
consideration further, and assuming that: a) all the clusters 
we have studied formed during the most recent starburst; and 
b) that during the starburst, 10% of stars formed in clus- 
ters that are still bound in the currently observable era, then 
we estimate the mass of the young field star population at 
~ 5 x 10 8 Mq. This value is an underestimate, as it was calcu- 
lated based on t he observable c lust er population. In addition, 
recent studies by lLarsenl (120091) and lGielesi (l2009h find that the 
underlying cluster mass function is likely to be a Schechter 
function, with a limiting mass of M» ~ 10 6 M Q . This is 
another factor that may lead to a false estimate of this total 
mass, as we have integrated over a power law with no im- 
posed maximum mass. While this value represents a lower 
limit, we note that the survival of a significant number of mas- 
sive clusters to ~ 150 Myr argues against a prolonged era of 
mass-independent cluster disruption. 

4.3. Age distribution 

Using the spectroscopically derived ages, we have plotted 
the age distribution of all clusters in the sample in Fig. |7J 
where the solid and dashed lines represent the disk and nu- 
clear region cluster populations respectively. We chose to bin 
the data in 50 Myr bins to reflect the errors in some of the 
measured ages; also, the first bin has half the size of the rest, 
to mirror the lack of disk clusters younger than ~ 40 Myr. 
We note that the overall shape of the histogram is maintained 
when a smaller bin size is applied, however the low number 
of clusters introduces uncertainties. The plotted error bars ac- 
count for Poissonian sampling and go as y/~N for each bin. As 
the spectroscopic data sample the entire galaxy disk and nu- 
cleus, this can be treated as a rough indicator of the M82 star 
formation history (SFH), where the number of clusters formed 
reflects the mean star formation rate (SFR) at that epoch. Nat- 



10 



I. S. Konstantopoulos et al. 



TABLE 3 

Spectroscopically derived ages and radial velocities for all clusters in the sample. 
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urally, the numbers cannot be treated as correctly 'to scale' as 
our dataset is not complete. 

The histogram confirms some of the 'milestones' we have 
discussed in our recent work on the M 82 cluster population 
(ISmith et alJ 120071: iBastian et all 120071 and Paper I). First, 
by placing the burst between 200 — 250 Myr, it reflect s the 
times cale for the last encounter with M8 1 at 220 Myr (lYunl 
119991) . the event that most likely triggered the disk starburst 
and stripped away part of the disk. Note that this marks the 
onset of the disk starburst, not the period of maximum for- 
mation that follows closely. This sets another observational 
constraint on th e age of the disk s tarburst (this was discussed 
in Paper I and ISmith et aT]|2007l, for region B). Finally, the 
first bin reflects the ongoing nuclear starburst. 

In addition, the age distribution peak for the full disk 
sample (that we place roughly at 150 Myr) is consistent 
w ith the 156 Myr di stribution peak for region B, as derived 
in lSmith et al.l d2007l) based on a sample of 35 UB VI age dated 
clusters. This supports our interpretation of region B as a reg- 



ular p art o f the galaxy disk (ISmith et al.ll2007l Paper I). 

In § !3.1l we discussed the possibility of a degeneracy in our 
cluster age dating technique. This is expressed through the 
existence of two best-fit troughs (occurring in the majority of 
cases), roughly placed at 200 ± 100 Myr and 600 ± 100 Myr. 
This degeneracy in line profile arises because of the evolution 
of the Balmer line profile both in strength and width. These 
line properties fluctuate in the first few hundred Myr, be- 
fore commencing a more or less monotonic evolution. More 
specifically, stellar Balmer lines gradually broaden until the 
age of ~ 300 Myr and become narrower once again at ~ 
600 Myr. This makes it difficult to distinguish between spec- 
tra either side of this 300 Myr mark using a statistical method. 
One deduction can therefore be made with certainty: that the 
clusters we observed are all younger than ~ 600 Myr. This 
adds to the findings of Paper I, where we concluded that the 
population of region M82-B is not as o ld as previously be- 
lieved (in excess of 1 Gyr, as presented in lde Griis et al.ll2003L 
and consequent works). This still does not exclude the possi- 
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FIG. 7. — Top: Distribution of cluster ages over the past 300 Myr (this 
range contains all observed clusters). This is a rough representation of the 
cluster formation history of M82, as the sample is incomplete. The solid 
and dashed lines traces cluster formation in the disk and nucleus respectively. 
As the spectroscopic sample spans the entire galaxy disk and nucleus, the 
number of clusters formed can be used as a diagnostic of the star formation 
rate, in a relative scale. The burst is placed at 220 Myr, and its peak SF epoch 
at ~ 150 Myr, after which the SFR gradually declines. After reaching an 
apparent post-starburst minimum at ~ 50 Myr, the disk starburst appears to 
halt, while the nuclear burst leads the SF in the galaxy. 
Bottom: Spatial distribution of ages for the spectroscopic cluster sample. 
The vertical lines delimit the nuclear region and the horizontal lines mark the 
median age for the region over which they span (i.e. western disk, eastern 
disk, nucleus, and entire sample in green). This shows no concentration of a 
certain age group in any region of the disk, indicating a global starburst event 
lasting the entire range of measured ages. 



bility of a broader age distribution for M82, possibly stretch- 
ing to ~ 600 Myr. 

In order to discern the better fit let us consider the exist- 
ing evidence. Fir st, the aforementioned study presented in 
I Smith et al. (2007) found the (independently derived) age dis- 
tribution to be extremely similar to the one presented in Fig. [71 
virtual inactivity prior to the burst at ~ 250 Myr and a grad- 
ual quenching in more recent times. In addition, this distribu- 
tion is centred at ~ 150 Myr and finds no clusters older than 
~ 300 Myr. Second, adopting the older age range in these 
cases results in an extended distribution with a gap between 
~ 300 — 500 Myr; this cannot be explained without employ- 
ing complex s cenaria such as multip le bursts 12 . In addition, 
recent work bv lMcOuinn et alj ( 120091) suggests that a starburst 
should last no longer than ~ 400 Myr, therefore advocating 
against an extended distribution for the starburst cluster popu- 
lation. Finally, a note on the method itself: the two-minimum 
structure does not always exist in both the Hf3 and H7 fits: 
For instance, cluster 1 12.2 (presented in Fig. 0) only shows a 
double minimum in the H/3 fit. In the H7 fit there is only one 
trough, describing the same age range as the H/3 fit. 

In light of this evidence, we adopt the younger ages for the 
all clusters where this degeneracy is present. 

Incompleteness is an issue that needs to be addressed in any 
study of star clusters in extra-galactic environments; an analy- 
sis that is applied to a fraction of a population can lead to mis- 

12 This, in turn, cannot be explained given the current observational evi- 
dence, as the burst was most probably caused by the last (and probably only) 
passage about M8 1 



interp retations (as treated in ISmith et ail 120071: iGieles et~aT] 
2007). In this case, the lack of clusters with age greater than 
~ 300 Myr could be the effect of such an observational bias. 
The completeness of our sample was estimated by inspecting 
the age-mass diagram of the clusters and by adopting a con- 
servative detection limit of my = 21 (see Fig. [6]). We find 
that if a significant population of clusters with ages between 
300 Myr and 1 Gyr and masses above ~ 10 4 M Q existed, it 
would have been included in our sample (assuming that they 
formed with the same intensity as their younger counterparts). 
We therefore confirm that the decrease in the number of clus- 
ters with ages greater than ~ 300 Myr (see Fig. |6| is real 
and not a result of incompleteness. In addition, the lack of 
older cluste r s in th e photometric (UBVI) sample presented in 
I Smith et at] d2007l where the limiting magnitude was fainter) 
also argues for a statistically viable sample. 

Thus we conclude that the full disk of M82 (mostly likely 
including the nucleus) entered an era of increased star/cluster 
formation rate no more than ~ 300 Myr ago, which contin- 
ued in the disk until at least r ~ 50 Myr, and is ongoing in 
the nucleus. Furthermore, no area in the galaxy, other than 
possibly the nuclear region, appears to have a cluster forma- 
tion history distinct to the rest of the disk. Unfortunately, we 
cannot derive the absolute cluster formation history for M82 
from thi s dataset due to sam ple selection and detection lim- 
its (e. g. iBastian et alJl2005h and uncertain kno wledge of the 
cluster disruption timesca le for this galaxy (e. g. IGieles et alj 
120051: lLamers et ai1l2005l) . The bottom panel of Fig. |7J com- 
bines the age and spatial distribution of clusters and clearly 
shows the scatter of all disk cluster ages about 150 Myr, and 
the very young ages of all nuclear clusters. It also shows the 
lack of young clusters at large galacto-centric distance: be- 
yond ~ 2 kpc, no clusters lie significantly below the distri- 
bution mean. This is unlikely to be a selection effect, as the 
brightest clusters were chosen for this programme, and these 
are most likely to be the youngest ones. 

4.4. On the star formation history ofM82 

We will now use the presented age measurements to inves- 
tigate the star formation history of M82. Although the sample 
is incomplete, star cluster ages can still be used as rough indi- 
cators of the CFR of the galaxy. 

The observational evidence for SF activity in the disk is 
as follows: there appears to be little or no recent (in the last 
Gyr or so) SF, prior to the onset of the starburst. The star- 
burst itself is observable in the increased cluster formation at 
t ru 200 Myr. The overall CFR in the disk appears to have 
dropped steadily since the peak starburst epoch at ~ 150 Myr, 
reaching an apparent hiatus ~ 50 Myr ago. However, given 
the errors inherent in this analysis (mainly due to incomplete- 
ness), we cannot firmly establish whether the CFR did indeed 
decrease or if it remained stable. Furthermore, our sample 
contains no evidence for the existence of young clusters in the 
outer regions of the galaxy disk, with younger disk clusters 
being situated closer to the nucleus. This deduction can be as- 
sumed to be true with some degree of confidence, as a young 
cluster forming region would shine bright in the outer regions, 
where the extinction is below average. Also, our observations 
were planned in such a way as to include the brightest (and 
therefore youngest) available sources; the ages of these clus- 
ters place their formation in the first ^100 Myr after the on- 
set of the burst. These observations pose questions about the 
mechanisms that maintained the M82 starburst over the past 
~ 250 Myr. 
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The nuclear burst appears at first glance to offer a simpler 
picture, with zero activity for most of the burst, and intense 
cluster formation in the current era. However, incompleteness 
effects are very grave in the nucleus, the result of both techni- 
cal/observational and physical factors: first, the nuclear region 
spreads over a mere ~ 200 pc, where only two slits could be 
accommodated. Second, older clusters will be outshone by 
the very youngest clusters and cluster-forming clumps. In ef- 
fect, the nucleus may have been producing clusters at a steady 
rate over the course of the disk starburst, as material was fun- 
nelled into the regions around the bar. Furthermore, the pre- 
cise location of the observed knots is not simple to pinpoint 
in this optical dataset. We discuss this further in the following 
section. 

Having considered this body of evidence, we propose the 
following scenario for the SFH of M82: following a pro- 
longed era of low activity in the disk, the SFR increased 
rapidly, probably triggered by the encounter with M81, some 
~ 220 Myr ago. As the starburst spread across the disk, 
the SFR increased and reached a maximum soon after (r ~ 
200 Myr). During this time (starting at the time of the burst), 
material star ted to slowly gravitate towar ds the nucleus (as 
proposed by iForster Schreiber et alJl2003l) . This withdrawal 
affected the outermost regions first and continued inwards, 
causing the SFR to drop gradually, as fewer and fewer re- 
gions were forming stars intensely. This process continued 
until t ~ 50 Myr, when disk activity dropped to a similar 
level as before the starburst. The nucleus may have entered 
an era of high SF along with the disk and maintained this high 
SFR throughout, to the present day, fed by the continuous sup- 
ply of material, channelled in by the bar. 

4.5. Cluster Radial velocities 

We present the measured cluster radial velocities (RVs; de- 
rived by cross-correlation with the best fitting spectroscopic 
models) and line-of-sight neutral interstellar gas RVs (derived 
from the Nal D doublet line at 5890/5895 A) in Fig. M For 
comparison, we have also plotted the major axis RV measure- 
ments derived from the near-IR Call stellar absorption and 
Pa(10 ) and [S ml nebular emission line s taken from the liter- 
ature dMcKeith et alJll993H Greve 2004). 

Our observations extend the major axis coverage of RV 
measurements significantly further than previous studies, al- 
though we caution that our data are not strictly 'major axis' 
measurements since our sample includes clusters located 
both above and below the disk mid-plane, at heights up 
to ~ 500 pc. The rotation curve, as traced by the Call, 
Pa(10) and [S III] data, is a composite of a steep inner sec- 
tion tracing the nuclear b ar dWills et al.l 120001; iGrevd 120041: 
IWestmoquette et ai]|2007l) . two flat (or gently rising) sections 
tracing the main disk (hereafter the flat part of the curve), and 
an intermediate transition region. 

In general, the cluster RVs are consistent with the Call 
measurements throughout the disk. This is in itself quite in- 
teresting considering the large range in heights above and be- 
low the disk sampled by our dataset. There are, of course, 
a number of notable exceptions: the off-disk cluster 152 
(at ~ —2000 pc) has a significantly lower RV than the stel- 
lar content and is discussed below as a candidate old popu- 
lous cluster. Clusters 108, 112.2, 125 and 126, located be- 
tween -500 and -900 pc, have RVs 60-100 kms -1 higher 
than the corresponding Ca II measurements, and neighbouring 
clusters. Besides their discrepant RVs, this group of clusters 
have no other distinguishing characteristics (their ages and 



colours are consistent with those of their neighbours), how- 
ever, they are associated with region B. Their high RVs there- 
fore support the hypothesis that region B represents a window 
of lower-than-average extinction, since if the clusters are lo- 
cated deeper within the galaxy, we would expect the radial 
component of their orbital velocity to be larger than that of 
material nearer the surface. We note that the discrepancy be- 
tween the Ca II and cluster RVs in region B implies that the 
Ca II absorption originates principally from the field stars (or 
their scattered light), not from the star clusters. 

Our observations include three slits in the nuclear re- 
gion, from which we have extracted six high S/N spectra. 
These probably represent bright 'cluster complexes' contain- 
ing blended light from a number of sources 13 . This is because 
our spatial resolution is not high enough to identify individual 
clusters within the crowded nuclear region. The RVs of five 
of these six cluster complexes place them on the steep part of 
the rotation curve, implying that they are associated with the 
nuclear bar. The exact location of these clusters with respect 
to the galaxy nucleus was briefly discussed in § 14.51 the evi- 
dence presented by the radial velocities suggest that they may 
not lie at the very deepest part of the galaxy nucleus - as also 
suggested by their relatively low extinction - therefore their 
characterisation as 'nuclear' may be inaccurate. Having said 
that, they provide the closest tracer to the behaviour and state 
of the nuclear region in our dataset. In order to study the true 
nuclear cluster population of this galaxy, IR observations are 
necessary. 

Interestingly, we find that, according to the Nal D veloc- 
ity measurements, the neutral gas does not appear to follow 
the same rotation pattern as the stars, clusters or ionized gas. 
Throughout the disk, the Na I D measurements are well de- 
scribed by a single, flat rotation curve with a gradient approx- 
imately equal to that of the stars/clusters in the outer disk. 
Nal D absorption must therefore originate from neutral gas 
located further out in the disk (given that they absorb clus- 
ter light), where the radial component of its orbital velocity is 
less. This also explains why none of the Nal D measurements 
are associated with the fast bar orbits in the nuclear region. 

4.6. Candidate Old Populous Clusters 

We now return to the two discussed OPC candidates. As 
discussed above, cluster 152 was found to have a lower-than- 
expected RV for its position, with a measurement more con- 
sistent with the neutral gas RVs, implying that it may lie at 
a considerable distance away from the disk. No certain RV 
measurement could be performed for the other OPC candi- 
date, no. 139, due to the low S/N of its spectrum (however a 
low-confidence fit places the cluster within the M82 system). 

Both clusters are located at projected heights below the disk 
of > 200 pc, implying that they may be part of the galaxy 
halo. A brief inspection of the area surrounding these clus- 
ters reveals the presence of numerous red stars, possibly red 
giants (RGs), or other, highly reddened stars. We investigated 
the possibility of photometric contamination by varying the 
aperture/sky annulus size and found no evidence of such an 
effect. 

Given the red colour and faintness of cluster 139 compared 
to the majority of clusters in the sample, there is a possibility 
that this object may in fact be a background galaxy - although 
we have no RV data to confirm or refute this. We can reject the 

13 As noted before, it is not possible to deduce whether these are physical 
associations or chance projections 
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FIG. 8. — Top panel: Cluster radial velocities (solid circles) plotted over the major axis rotation curve de fined by the near-IR C a II s tellar absorpti on (upward- 
facing triangles) and [S III] and Pa(10) nebular emission measurements (downward-facing triangles) from McKeith et al. ( 1993) and Greve (2004). The cross 
indicates the position of the dynamical centre of M82. Rectangles represent the radial velocity of the Na I D doublet absorption line (tracing neutral interstellar 
gas) measured in each cluster spectrum. As these measurements do not require a high S/N, we have obtained a ti^a I D for every slit in the sample, leading 
to Nal D data with no associated cluster velocities (albeit with larger errors). Bottom panel: the reddening distribution of the spectroscopic cluster sample, 
as derived using broadband colours, and corrected for foreground emission from the MW. The dashed vertical lines roughly outline the nuclear region and the 
horizontal lines indicate the median extinction value for each region. The blue boxes represent the OPC candidates. 



possibility of cluster 152 being a background galaxy (regard- 
less of its photometry), based simply on its radial velocity, 
which places it within the M82 system. In § |3.1| we found evi- 
dence that this cluster is older than ~ 8 Gyr (using a |Z fit). 
Given the uncertainty of the Padova isochrones for interme- 
diate/old ages, we assume that this cluster is an old Globular. 
However, the possibility still exists that cluster 152 is of in- 
termediate age, constituting the first potential detection of the 
sort in M82. In any case, this object deems further study. 

5. PROJECTION EFFECTS: A PROPOSED TOY MODEL FOR THE 
INVISIBLE PHYSICAL STRUCTURE OF M82 

A careful study of all the results we have presented in this 
paper can reveal many pieces of information on the physi- 
cal structure of the galaxy, which is obscured in our nearly 
edge-on viewpoint. Such results have been presented in the 
past, but they have never been at the focus of any study. It 
has been discussed that M82 is a spiral galaxy (as far back 
as OM7 8), with a nuclear bar spanning some 170 pc in pro- 
jection (lAchtermann & Lac vl [19951: IWills et al.ll2000l) . which 
extends along our line of sight. Observations in Ha and in 
the IR have also allowed for a detailed study of the super- 
galactic wind emanating from the galaxy nucleus, fuelled by 



the nuclear clusters dWestmoquette et alj|2007ll2009l) . In this 
section we will present a model for the structure of the galaxy, 
based on the mentioned observations and the ones presented 
in this paper. 

5.1. Modelling disk kinematics and matching observed 
velocities 

Our results find a distribution of radial velocities typical of 
a disk galaxy, with solid body rotation in the nucleus (slits 67 
and 78) and flat (or gently outwards increasing) components 
for disk velocities. There are, however, a handful of clusters in 
region M82-B (200 to 1000 pc east of the nucleus) whose ve- 
locities do not follow the observed trend: clusters 108, 1 12.2, 
125 and 126 all have significantly higher velocities than their 
neighbours and the underlying stellar component. A simple 
consideration of M82 as it would be seen face on reveals that 
clusters situated at similar projected galactocentric distance 
may reside at a range of depths along the line of sight. A rep- 
resentation of this effect may be found in Fig. [9] left panel, 
where we have assumed that the galaxy disk is roughly circu- 
lar. In this schematic, the thick line represents the nuclear bar, 
and crosses map out a grid limited to a 3 kpc radius (the pro- 
jected visible radius of M82); vectors represent the expected 
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radial velocity of an object in differential rotation at that lo- 
cation. We note here that this is a simplified representation of 
the galaxy, where the intricate dynamics of the regions near 
the bar are not modelled, and only regular, circular disk ro- 
tation is considered. This is simply because there is very lit- 
tle information to be revealed by nuclear clusters concerning 
galaxy kinematics. 

This basic trigonometric consideration reveals that a clus- 
ter on a regular orbit near the visible surface of the galaxy 
will appear to move significantly slower than another regu- 
larly orbiting cluster situated at great depth, r. The weight of 
this effect ranges with projected galactocentric distance (in- 
dicated by d in Fig. at the extremes of the galaxy disk 
this difference will be within the measurement errors; in the 
mid-regions of the galaxy, assuming the M82 disk is roughly 
circular, the velocities will differ by a factor of two. Clusters 
projected on top of the nucleus will appear to have near-zero 
RVs, as their velocity is mostly in the tangential direction and 
is therefore not measureable; just to the side of the nucleus 
the discrepancy reaches a factor of 10. 

5.2. Spiral structure 

Any existing spiral structure in M82 will affect its appear- 
ance, even in our highly inclined viewpoint. Given the orien- 
tation of the bar, which extends along our line of sight, one 
would expect spiral arms to originate from the visible central 
regions of the galaxy on the near side, and a corresponding 
site in the obscured side. This could result in a highly irregu- 
lar dust distribution, and therefore significant discrepancies in 
the visible components of the disk. 

In order to find evidence of the two-armed, b arred spi- 
ral st ructure (as hypothesised in the literature, see lYun et all 
11994 and references therein) one must also consider the in- 
clination of the disk, which provides a view of M82 from be- 
low. Assuming no significant spiral sub-structure other than 
the two dominant arms, we should expect to face one strong 
extinction front in the one side of the galaxy (the near spi- 
ral arm), and an 'opening' in the dust distribution on the 
other side, along lines of sight underneath the bulk of the 
galaxy, and into the disk. More specifically, the near spiral 
arm should stem f rom the tip of the bar (the approximate lo- 
cation of SSC Al I Smith et al. 2006) and extend westwards, 
and the far-side arm should originate at a symmetric location 
an d 'wrap' around reg ion B. Such a conclusion was reached 
by lMayya et al.1 (120051) . who derived the approximate location 
of the spiral arms by subtracting a hypothetical axisymmetric 
exponential disk from NIR imaging of the galaxy. 

This structure is observed in the optical and our hypothe- 
sis on the approximate placing of the two spiral arms is sup- 
ported by long wavelength observations (IR and radio). Ex- 
amining the optical data presented in this paper, we assign the 
proposed structure as follows: the western side of the galaxy 
hosts the near spiral arm, hence the highly obscuring extinc- 
tion front in the southwest disk. A faint trail of reddening 
dust can be seen superimposed on the central region, stem- 
ming from the general vicinity of tip of the nuclear bar, and 
spreading towards this extinction front. Therefore, the eastern 
side is expected to host the far spiral arm; as it hosts region B, 
where kinematics reveal lines of sight deep into the body of 
the galaxy, it is consistent with this concept. 

Longer wavelength studies lend additional support to 
this model. First, we co nsider Hubble NICMOS data 
dAlonso-Herrero et al.l 1200 lb in the J and H near-infrared 
bands, which trace star light, with much less obscuration from 



interstellar dust . These data show the eastern disk as a 
crowded stellar field, with a large number of star clusters that 
are invisible in sub-micron wavebands. The part of region B 
closest to the nucleus (200 to 700 pc) is more crowded, as 
it represents a cross section of more than ~ 4 kpc in depth, 
and projections of stars and clusters are the norm. Unfor- 
tunately, no equivalent high-resolution IR data exist for the 
western disk. 

As a tracer of spiral arms we have sought out radio data with 
continuou s coverage of the g alaxy. We use the CO molecu- 
lar data in W alter et al.l d2002l see their Fig. 2 where radio and 
optical data are overlayed) as an appropriate diagnostic of spi- 
ral structure - CO is a tracer of H2. We find a good agreement 
between the structures we propose as repr esenting the spiral 
arms and the CO streamers discovered by I Walter et al.l The 
radial velocities found for these streamers are consistent with 
material delineating the outer parts of the disk, as a spiral arm 
would. 

This model finds support in our analysis of disk kinematics; 
in the right panel of Fig. [9] we have matched observed clus- 
ter velocities onto the modelled RV field, according to their 
measured projected galactocentric distance (d). Applying this 
method, we find clusters in region B to be consistent with be- 
ing located at large radii (r), given their high radial velocities. 

Finally, we consider the placement of LOS gas, traced by 
the Na I D doublet line (plotted onto Fig.[9]as filled red boxes). 
We expect Na I D line velocities to have maxima at extreme d 
(i. e. the western and eastern edges of the disk), and a linear 
trend between those extrema crossing (d, r) = (0, 0); this is 
confirmed by our observations. There is one deviation from 
this, as we find the interstellar gas to be situated at greater 
depth across region B (but at the surface between r = 600 — 
800), thus lending some more support to our hypothesis. 

Na I D measurements can also provide more insight on the 
exact location of their corresponding clusters. As it is not pos- 
sible to distinguish kinematically between symmetric near/far 
side locations, we resort to indirect means to choose between 
the two solutions to this problem: in some cases, the Nal D 
clump is placed behind the respective cluster (wrt our LOS), 
which is impossible as the line occurs in absorption of clus- 
ter light. In these cases we deduce a far side location for the 
cluster. 

We therefore present new observational evidence for the 
morphological structure of M82, and confirm that it consists 
of two dominant spiral arms a nd very weak, if any , other 
spiral structure, as deduced by llchikawa et al.l (1 19951) using 
NIR im aging of the galaxy, and followed up by iMayya et al.l 
(120051) . 

5.3. Implications of the proposed model on region M82-B 

Under the proposed model, part of the eastern M82 disk is 
left relatively free of extinction. The spiral arm is situated at 
the far reaches of the galaxy, therefore all owing a clear view 
into the body of the galaxy. In Paper I and lSmith et al.l (120071) 
we argued that region B in the eastern disk is seen through 
low extinction windows, similar to "Baade's windows" in the 
Milky Way. Here, we offer the interpretation that these win- 
dows form as a result of the inclination of M82 and provide 
a line of sight under the spiral arm extinction front. Thus, we 
confirm that region B is not intrinsically brighter, but simply 

14 Central wavelengths are 1.10 and 1.60^tm for the J and H bands re- 
spectively; assuming a galactic extinction law, the exact relation is Aj = 0.3 
and An = 0.1 mag for each magnitude of extinction in V. 
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FIG. 9. — Left: A model of the M82 RV-field. The orientation is east to the left, and north out of the page (following our observed viewpoint of M82). Crosses 
represent grid positions with galactocentric distance less than 3 kpc, which we assume to be the full extent of the (hypothetically circular) M82 disk. The thick 
line down the centre of the plot represents the approximate location of the bar (the extremes on the r-axis are uncertain, as we can only observe its projected 
location). Vectors represent modelled radial velocities, and are colour-coded as receding (red, east side) and approaching (blue, west side). The grid locations 
with plotted velocities are all outside the inner, solid-body rotation region associated with the bar (represented by the thick tilted dark line). The observable radial 
component of this velocity ranges from zero at d = 0, where the velocity is tangential to our line of sight, to the maximum velocity on the r = line (where all 
of the motion is in the radial direction). Right: The same model, with overplotted predicted cluster locations (stars) based on the kinematics of the cluster and 
nearest pocket of interstellar gas (Na I D, red squares). This is not an accurate depiction, rather a diagnostic to confirm our hypothesis on the location of clusters 
in region B, where this demonstrates the depth at which these clusters are located. 



seen under a clearer, relatively unextinguished line of sight. 
This interpretation is supported by the high radial velocities 
of the mentioned group of region B clusters, as they lie within 
these low extinction windows. This means that they can be 
seen clearly, even though they are situated deep within the 
body of the galaxy, and they do not form a distinct population, 
in terms of SFR, intrinsic luminosity, or actual kinematics. 

6. SUMMARY 

In this paper, we have presented a spectroscopic and pho- 
tometric study of the disk cluster population of M82, using 
a sample of 44 bright, isolated disk clusters. We have also 
studied five cluster complexes in the galaxy nucleus, which, 
together with the 44 disk clusters, provide virtually contin- 
uous coverage across the galaxy and comprise the largest to 
date spectroscopic sample of extragalactic young clusters. 

Our primary goal was to derive spectroscopic (age, radial 
velocity) and photometric (colour, extinction, mass) informa- 
tion on the clusters, and to characterise the population based 
on the properties of the individuals. In addition, we have ex- 
tracted some information on the environment which they in- 
habit, in terms of the M82 starburst history. 

Having combined information from spectroscopy and 
imaging, we found the disk clusters to form a uniform pop- 
ulation, displaying no positional dependencies with respect to 
age. We have therefore roughly derived the starburst history 
of M82, based on the cluster age distribution: this we describe 
as a long period of low SF activity, followed by a burst of star 
formation that gradually decreased in the disk, while the nu- 
cleus has been starbursting for at least the last 15 Myr, and 
perhaps throughout the starburst, and appears to be leading 
SF in the galaxy at the current epoch. 

In summary, we conclude the following: 

• The disk cluster population of M82 is consistent with 
having formed during the last encounter with M81, 
~ 220 Myr ago. This is based on spectroscopic ages, 



which range between (30 — 270) Myr. 

• The SFR of M82 reached a peak some ~ 150 Myr 
ago, and has remained virtually constant since, with 
the disk starburst leading star/cluster formation up to 
t ~ 50 Myr, and the nuclear starburst taking over since. 

• The nuclear starburst appears to be ongoing, as indi- 
cated by the possible detection of WR features in the 
central regions of the galaxy. 

• The only possible evidence of pre-starburst cluster for- 
mation comes from cluster 152, a candidate old popu- 
lous / Globular cluster, that was observed for the first 
time. 

• Cluster reddening ranges between 0.4 — 3.2 mag in the 
disk and 1.4 — 3.5 mag in the nuclear region. The mea- 
sured masses range between typical low-end values of 
~ 10 4 M and reach ~ 10 6 M Q . From these values, 
a total mass of ~ 10 8 is deduced for the star cluster 
population of M82. 

• The cluster system follows a regular rotation pattern 
about the nucleus, with a nearly flat rotation curve com- 
ponent for the disk, and a solid-body rotation part in the 
bar area. 

• Region M82-B is a regular part of the galaxy disk, and 
represents a line of sight into the body of the galaxy 
through lower average extinction windows, made possi- 
ble by the inclination and spiral structure of the galaxy. 

• The data presented in this work support previous pre- 
dictions of a two spiral arm system, as indicated by our 
study of cluster radial velocities. We place the roots of 
the near spiral arm in the vicinity of SSC Al (the tip of 
the bar) extending west, and the end of the far-side arm 
above region B in our subtended viewpoint. 
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This concludes the spectroscopic and photo metric study 
of the M82 cluster population we initiated in I Smith et all 
J20071) . This series has included, apart from Paper I, in- 
vestigations of the enviro nment of the M82 starburst core 
dWestmoqu ette et al. 200 7]) and the dy namics of the disk and 
inner wind dWestmoquette et aLll2009l) . 
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